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^ ' Abstract 

A strangeness and isospin asymmetry including statistical bootstrap model anal- 



ysis of the multiparticle system produced in the Pb + Pb collision at 158 AGeV at 
CERN is presented. It is concluded that this interaction process has not crossed the 
deconfinenent line. Direct comparisons with the results of similar analyses pertaining 
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to nucleus-nucleus and p — p collisions at CERN are made. The overall picture points 
to the S + S collision at 200 AGeV as a prime candidate of a process which has crossed 
the border separating the hadronic from the deconfined phase of matter. 
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The quintessential issue probed by relativistic heavy ion experiments concerns whether 
or not a phase of matter beyond the hadronic one has been achieved during the interaction 
process. Any conclusive evidence, furnishing an affirmative answer to this question, is rel- 
evant not only in solidifying the status of QCD, as the microscopic theory for the strong 
interaction, but also in providing a short glimpse into the early cosmos. 

The most tangible element of experimental input, on which one relies to infer as to 
what has taken place during a given collision between two heavy nuclei, is the constitution 
of the produced multiparticle system. One viewpoint to examine is whether such a system 
originates from a thermally and chemically equilibrated state, having gone through a complex 
stage of dynamical evolution. This is the premise upon which a number of extensive analyses 
have been conducted [1-3] , pertaining to multihadron systems in heavy ion collisions, as well 
as p — p and e + — e~. A common feature in the aforementioned studies is the adoption of an 
"ideal hadron gas" model to simulate such multiparticle systems. These kinds of approaches, 
however, can go only so far as to verify the thermalization hypothesis - which, in fact, they 
do - but fall short of addressing the basic question whether a given observed multiparticle 
system provides any signature which could trace its source inside or outside the hadronic 
phase of matter. This is due to the absence of interactions in the thermodynamic description 
of such relativistic, many-particle system in this type of calculations^. 

The fundamental issue of constructing a self-consistent scheme, which furnishes a ther- 
modynamic account of an interacting relativistic multihadron system was formulated in the 
1960's by Hagedorn [4-7] and is known as "Statistical Bootstrap Model" (SBM). In this 
scheme, the notion of interaction at a distance between relativistic particles is simulated 
by a successive organization of the system in terms of particle-like entities of increasing 
complexity, known as fireballs. At the lowest level enter all known hadrons, which are in- 
corporated as input into a self-consistent (bootstrap) iterative scheme, realized in terms of 
an integral equation for the fireball mass spectrum function p(m). The final construction of 
the SBM is accomplished by enriching the scheme, in a relativistically consistent way, via 
the introduction of thermodynamic variables (temperature and fugacities), which monitor 
its thermal and chemical equilibrium. After some interim mathematical manipulations, the 
1 Only an excluded volume correction has been addressed in some of these studies. 
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bootstrap equation acquires the form 

<p(T, {A}) = 2G(T, {A}) - exp[G(T, {A})] + 1 , (1) 

where f is the so-called input function, whose specification is determined by masses and 
quantum numbers of the known hadrons, while G incorporates the fireball mass spectrum. 
The most notable feature of the bootstrap equation is that it displays, in the f — G plane, 
a square-root branch point at 

V?(T cr , {A cr }) = In 4 - 1 , G(T cr , {A cr }) = In 2 . (2) 

The physical branch serves to define a critical hyper-surface in the space of thermodynamic 
variables, which sets the limits of the hadronic phase. Finally, the "statistical" aspect of 
the SBM is encoded in a grand-partition function Z$bm- It is of crucial importance to note 
that the mass spectrum function, which incorporates the interaction dynamics among the 
hadrons, enters the expression for Z$bm- 

The task of employing the SBM to perform non-trivial (in the sense of including in- 
teractions) thermal-type analysis of multiparticle yields from heavy ion collisions has been 
undertaken by the present authors in recent publications [8-11]. To this end, we extended 
the original Hagedorn scheme, in which a single fugacity variable for baryon number was 
employed, so as to include fugacities relating to strangeness numbers, both net (strangeness 
minus anti-strangeness) A s and absolute (strangeness plus anti-strangeness) A|s| [9-10], as 
well as a "net charge" fugacity [11]. The latter, in combination with the baryon number 
fugacity, accounts for isospin asymmetry in the colliding system. We formulated in this way 
the Strangeness-including SBM (SSBM). Moreover, we departed [8-9] from the usual choice 
of a certain parameter a entering the bootstrap scheme, which determines the particular 
partitioning between a kinematical and a dynamical factor composing the mass spectrum 
function G. In this way, a direct relationship arises between the maximum value of the tem- 
perature, To, on the critical surface and the MIT bag constant. The SSBM is analytically 
applicable only within and up to the limit of the hadronic phase, defining in a precise way 
this limit. For our purposes the boundary of the hadronic domain will be designated by its 
projection on the 2-dimensional {T,ji u ) phase diagramme. This projection is close to the 
intersection of the critical surface with the ji s = plane. 
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Given the wide interest generated by the Pb + Pb experiments at SPS, we shall present, 
in this letter, results stemming from an SSBM analysis of data at 158 AGeV. We shall also 
make a global assessment of all results we have so far obtained via the SSBM approach, 
which include S + S, S + Ag and p + p collisions at SPS. Such an analysis is provided by 
the space of thermodynamical variables, a comprehensive form of which is represented by 
the set (V, T, A M , A^, A s , j s ), where 7 S = A|s|. Two constraints are imposed on these variables, 
corresponding to zero strangeness and the connection between total baryon number and total 
charge of the colliding system 

<Q > N™ w 

where N™ (N^ 1 ) is the total number of participant protons (neutrons). 

Referring to the first column of Table 1, we consider the full phase space yields of particles 
produced in the Pb+Pb NA49-experiment at CERN [12,13] and perform a x 2_ fit i n relation to 
the theoretical expressions furnishing the corresponding particle multiplicities in the SSBM 



Nfh =x din Zssbm 



(4) 

{Ai=l} 



through which an optimum value for each one of the thermodynamical variables is deter- 
mined^]. 

Columns 3 and 4 of Table 1 present the SSBM calculated particle yields corresponding to 
optimum values of the thermodynamical variables, obtained with and without the inclusion 
of pions, respectively. The corresponding optimum values of these variables are entered 
in Table 2. Finally, in Table 3 we list the experimentally measured particle ratios, upon 
which we shall base the presentation of our results. In Fig. 1 we depict the location of the 
multiparticle emitting source on the T — \x u plane, resulting from the SSBM analysis with 
the inclusion of pion multiplicities and the fugacity variable 7 S fixed to its optimal value 0.76 
(see second column of Table 2). For each particle ratio we have formed corresponding bands, 
as determined by experimental uncertainties. The source coordinates, obtained from the fit 

to the particle yields, are located at the point which is nearer to all the particle-ratio bands, 
2 Actually, along with the thermodynamical variables, two additional ones enter the x 2_n t, namely the 
two Lagrange multipliers that serve to impose the constraints. 
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identifying the region of the original achievement of thermal and chemical equilibrium by 
the multiparticle system. The most significant feature in the plot is the thick solid line, the 
SSBM (de) confinement line for the value of 7 S corresponding to the Pb + Pb interaction. 
It is an indigenous characteristic of the SSBM and its particular specification in the plot 
involves the choice for the critical temperature, T , at zero chemical potentials /i^ and /i s 
(7 S is fixed). The reasoning behind our particular choice: To — 183 MeV, corresponding to 
maximum MIT bag constant value B 1 ^ = 235 MeV [12] is that it allows as much as possible 
space to the hadronic phase [9]Q Thus no doubt can be left that an interaction point is 
outside the hadronic domain if it lies beyond the surface for which the hadronic domain is 
maximally extended. 

Let us point out that the recent data [14] give smaller values for the An multiplicities of 
S + and for the A/A ratio. This leads to the lowering of the fitted temperature as well as 
the lowering of j s . If we had used the older values H~ + S + = 8.19 ± 1.06, H~ = 7.23 ± 0.88 
[15] and A/A = 0.20 ± 0.04 [16] we would get T = 172.3 ± 9.3 MeV, fx u = 70.4 ± 7.4 MeV, 
H d = 80.3 ± 7.8 MeV, /i s = 10.4 ± 9.7 MeV and 7, = 0.827 ± 0.088 with x 2 /dof = 22.8/7 
and even larger values for T and 7^ if the pions were excluded. These results are compatible 
with the results of Analysis B for the Pb + Pb case in [16]. 

From Fig. 1 it is evident that, the source of the produced multiparticle system in the 
Pb + Pb interaction at 158 AGeV lies well within the hadronic domain, which equivalently 
means that the deconfined, partonic phase has not been attained in this system. 

It is evident when comparing the last columns of Table 2 that neither a notable improve- 
ment of the fit (e.g. the value of x 2 /dof) nor a remarkable change in the fitted thermody- 
namic variables occurs when the < it > multiplicity (which contains pions) is excluded from 
the data, as is the case with the fits in the S + S [10] and S + Ag [11] data. The same 
conclusion can be inferred from Fig. 2, where the fittings of particle multiplicities in 4n 
phase space, with and without the inclusion of the pion multiplicities, are compared to the 
experimental values. The overall picture is that there is no notable difference between the 
two cases, meaning that the produced entropy (mainly associated with the pion production) 
3 Also To(/ig = 0) ~ 183 MeV is the maximum temperature for non- negative /i s in the HG domain. Recent 
lattice QCD calculations give Tq ~ 175 MeV for [2+1] quark flavours [13]. 
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is well accounted for by the hadronic thermal model. 

A summary of results attained through the SSBM analysis of S + S [10] and S + Ag 
[11] interactions at 200 AGeV (NA35), as well as the p + p (UA5) at several energies [10], 
is exhibited in Fig. 3. It should be pointed out that the surface setting boundaries on the 
hadronic domain is a constraint among the thermodynamic variables (T, A u , A<j, A S ,7 S ). The 
projection of this surface on the (T, /ig) plane which corresponds to the particular value of 7 S 
resulting from the fit to the experimental multiplicities is plotted for every interaction. Thus 
every point of interaction should be compared with the relevant projection. The analysis 
shows that the S + S source is situated mostly outside the hadronic phase (with probability 
of being outside 74%), whilst the S + Ag is located just beyond the deconfinement line (with 
probability of being outside 52%)n. In addition, the analysis exhibits a large (~ 30%) entropy 
enhancement of the experimental negative-hadron yields compared to the model, an effect 
observed also by other calculations [1-3]. This enhancement may be attributed to contribu- 
tions from the deconfined quark phase with many libarated new partonic degrees of freedom. 
Also, the thermal-statistical models [1-3], which do not posses inherently any bounds for the 
hadronic phase, give fitted source temperature beyond those of the corresponding SSBM 
analysis: T ~ 180 — 200 MeV, which is about 5-25 MeV higher. 

We conclude that all these results are corroborating our position that the Pb + Pb in- 
teraction at the maximum energy of 158 AGeV is well within the hadronic phase. The 
S'-induced interactions and in particular the S + S one are located beyond the HG phase, 
in the deconfined quark-gluon domain. For the p + p collision at ^fs = 200 to 900 GeV, the 
SSBM analysis clearly locates this system within the hadronic domainQ 

In summary, the potential of SSBM to perform not only thermal-statistical fits to mul- 

tiparticle systems, but also to define the limits of the hadronic phase has been employed to 

It should be noted that, as far as experiments whose multiparticle data flirt with the critical curve, 

special considerations of computational nature are called for. Such issues have been discussed at length in 

Refs [10,11]. Moreover, as an additional aid to probing the region beyond the critical surface we conducted 

relevant analyses for an increased value of T . 

5 The p + p thermodynamic variables have been obtained through a grand-canonical analysis. It is known 

that canonical suppression is relevant to the p + p interaction and it may affect the extraction of parameters. 

It cannot, however, move the point outside the critical surface since it lies well inside the hadronic domain. 
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analyze and asses several nucleus-nucleus interactions. On the basis of these analyses we 
reach the conclusion that the Pb + Pb interaction at 158 AGeV has not crossed the decon- 
finement line, remaining within the hadronic phase. On the other hand, the S + S interaction 
at 200 AGeV appears to have produced a deconfmed partonic state, for the first time. 
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Table Captions 

Table 1 Experimentally measured full phase space multiplicities in the NA49 Pb + Pb ex- 
periment at 158 AGeV and their theoretically fitted values by SSBM, with the inclusion of 
the < 7r > multiplicity and without it. 

Table 2 Results of the analysis by SSBM of the experimental data from Pb + Pb experiment 
(47T phase space), with the inclusion of the < n > multiplicity and without it. 

Table 3 Particle ratios from the experimentally measured full phase space multiplicities in 
the Pb + Pb experiment at 158 AGeV, used in the analysis. 

Pb + Pb (NA49) Full phase space 



Particles 


Experimental 


Calculated 


Calculated 




Data 


With < 7T > 


Without < 7T > 


< 7T > a 


600 ± 30 


570.51 6 


527.05 c 


K + 


95 ± 10 


96.583 


93.470 


K~ 


50 ±5 


58.168 


54.528 


K s ° 


60 ± 12 


76.110 


72.800 


V 


140 ± 12 


150.73 


148.01 


V 


10 ± 1.7 


8.1956 


7.2435 




7.6 ± 1.1 


7.4101 


7.3914 




4.42 ±0.31 


4.0838 


4.1886 




0.74 ±0.04 


0.76299 


0.77076 


B-B 


362 ± 12 


366.65 


364.43 


A 


5.14 ±0.6 


4.8591 


4.6047 


A 


52 ±3 


48.263 


48.472 



a < 7T >= (7T+ +0/2 

b A correction factor 1.03272 has been included for the effect of Bose statistics. 

c A correction factor 1.03107 has been included for the effect of Bose statistics. This multiplicity is not 
included in the fit. 
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Table 1 



Pb + Pb (NA49) Full phase space 



Fitted Parameters 


Fitted with < tx > 


Fitted without < n > 


T (MeVl 

J. 1 J.VXC V I 


1 ^6 3 ± 4 2 


1 57 4 + 4 3 

XtJ I . t: I t: . -J 


\ 

11 


1 606 ± 051 


l.UtJtJ _l_ U.uUt: 


A .J 


1 686 ± 05Q 


1 721 ± 077 


\ 

^s 


1 1 78 ± 030 


1 171 +0 031 


7s 




U. I C/U _L_ U.UOX 


VT 3 /4n 3 


15.7 ±2.8 


13.8 ±3.4 


X 2 /dof 


10.85 / 8 


8.46 / 7 


Hu (MeV) 


74.0 ± 5.3 


77.2 ±6.6 


^ (MeV) 


81.7 ±5.9 


85.5 ±7.5 


/is (MeV) 


25.6 ±4.0 


24.9 ±4.2 



Table 2 



Pb + Pb (NA49) Full phase space 



Particle Ratios used 


Experimental 


for the plots 


Values 


< 7T > /{B - B) 


1.66 ±0.10 


K+/(B-B) 


0.262 ±0.029 


K~/(B-B) 


0.138 ±0.015 


K S °/(B-B) 


0.166 ±0.034 


p/(B-B) 


0.387 ±0.036 


p/(B-B) 


0.0276 ±0.0048 


<f>/(B - B) 


0.0210 ±0.0031 


E-/(B-B) 


0.0122 ±0.0009 


E + /(B-B) 


0.00204 ±0.00013 


A/(B-B) 


0.0142 ±0.0017 


A/(B-B) 


0.144 ±0.010 



Table 3 



Figure Captions 

Figure 1 Experimental particle ratios in the (/i n , T)-plane for the Pb+Pb interaction mea- 
sured in 4ir phase space with 7 S set to 0.76. The point and the cross correspond to the x 2 
fit with the < 7r >. The thick solid line represent the limits of the hadronic phase (HG) as 
set by the SSBM. The smallest experimental value of the ratio p/(B — B) — 0.351 cannot 
be depicted on the (/i u — T) plane because it has no solution for the given variables. Thus 
the space which is compatible with the experimental values is the one which is enclosed by 
the ratio with the largest experimental value. 

Figure 2 Comparison between the experimentally measured multiplicities in 47r phase space 
and the theoretically calculated values in the fit with < 7r > and without < ir > for the 
Pb + Pb interaction. The difference is measured in units of the relevant experimental error. 

Figure 3 (//#, T)-phase diagramme with points obtained from fits to p + p, S + S, S + Ag 
and Pb + Pb data and corresponding critical curves given by SSBM. 
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Pb + Pb , full phase space, y s =0.76 
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